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Abstract 

Recently Hofava proposed a renormalizable gravity theory with higher spatial derivatives in four 
dimensions which reduces to Einstein gravity with a non-vanishing cosmological constant in IR but 
with improved UV behaviors. Here, I consider a non-trivial test of the new gravity theory in FRW 
universe by considering an IR modification which breaks "softly" the detailed balance condition in 
the original Hofava model. I separate the dark energy parts from the usual Einstein gravity parts 
• in the Friedman equations and obtain the formula of the equations of state parameter. The IR 

Q . modified Hofava gravity seems to be consistent with the current observational data but we need 

\ some more refined data sets to see whether the theory is really consistent with our universe. From 

the consistency of our theory, I obtain some constraints on the allowed values of wq and w a in the 
Chevallier, Polarski, and Linder's parametrization and this may be tested in the near future, by 
sharpening the data sets. 
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Recently Hofava proposed a renormalizable gravity theory with higher spatial derivatives 
(up to sixth order) in four dimensions which reduces to Einstein gravity with a non-vanishing 
cosmological constant in IR but with improved UV behaviors by abandoning Einstein's 
equal-footing treatment of space and time H, @]. Since then various aspects have been stud- 
ied. In particular, in [jj], it has been pointed out that the black hole solution in the Hofava 
model does not recover the usual Schwarzschild-AdS black hole even though the general 
relativity is recovered in IR at the action level. (For another problem in cosmology, see 
For this reason, in [EJ an IR modification which allows the flat Minkowski vacuum has been 
studied by introducing a term proportional to the Ricci scalar of the three-geometry fx 4 R^ 
(for related discussions, see also [6J) and recently the general black hole and cosmological 
solutions have been found 0] , which reduce to those of [3j in the absence of the IR modifica- 
tion term and those of [E| for vanishing cosmological constant (~ A^y). (For other aspects, 
see 0). 

On the other hand, in pjj, the author argued that the dark energy may be explained by the 
Hofava gravity and obtained the equation of state parameter which seems to be consistent 
with the observational data, by neglecting the matter contributions 1 . 

In this paper, I consider an improved analysis of the proposal by comparing with the 
latest data which does not need to know about matter contributions, separately. This would 
provide a possible test of Hofava gravity. 

To this ends, I start by considering the ADM decomposition of the metric 

ds 2 = -N 2 c 2 dt 2 + g i3 (dx l + N l dt) (dx ] + N j dt) (1) 
and the IR-modified Hofava gravity action which reads 



S g = J dtd^N ± faK* - \K 2 ) - ^C^C* + ^ k B^V jR m k 

8 10 8(3A - 1) I 4 v ; w w 8(3A - 1) 



where 



Kij = ^ (9ij - ViNj - V d Ni) (3) 
is the extrinsic curvature (the dot (') denotes the derivative with respect to t), 

C ij = e iM Vk ^ R (3) j{ _ 1^(3)5^ ( 4 ) 

is the Cotton tensor, k, A, z/, /i, A^, uo are constant parameters. The last term, which has 
been introduced in 0,0,0], represents a "soft" violation of the "detailed balance" condition 
in [2J1 and this modifies the IR behaviors. 



Recently, it has been also proposed by Mukohyama l9| that the dark "matter" as integration constant in 
Hofava gravity, from the "projectability" condition [2|. But in this paper, I consider the non-projectable 
case, for the definiteness. 
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Now, in order to study the cosmological implications of the action (j2J), I consider a 
homogeneous and isotropic cosmological solution with the standard FRW form 2 



ds 5 



-c 2 dt 2 + a 2 (t) 



dr 2 



kr 2 /Rl 



(de 



2 , -2 

+ sin 



(5) 



where k = +1, 0, —1 correspond to a closed, flat, and open universe, respectively, and R is 
the radius of spatial curvature of the universe in the current epoch. Assuming the matter 
contribution to be of the form of a perfect fluid with the energy density p and pressure p, I 
find that Q 



K 



6(3A- 1) 

K 2 

6(3A- 1) 
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3/^/i 



8(3A- r 
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2k{A 
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--(p + 3p)± 



/ —k 2 

3^V 
8(3A- 1) \R$a 4 



R 2 a 2 
k 2 



(6) 
(7) 



where the prime (') denotes the derivative with respect to r. I have considered the analytic 



continuation [i A 



-fi 2 for Aw > (1, 0, EH and the upper (lower) sign denotes the 
Aw < (Aw > 0) case. Note that the 1/a 4 term, which is the contribution from the 
higher-derivative terms in the action (T5]), exists only for k ^ and become dominant for 
small a, implying that the cosmological solutions of general relativity are recovered at large 
scales. The first Friedman equation ([6]) generalizes those of jjjf and fBI to the case with an 
arbitrary cosmological constant and the soft IR modification term in [2, 13, E|- However, it is 
interesting to note that there is no contribution from the soft IR modification to the second 
Friedman equation ([7]) and this is identical to that of [3j. 

If the Friedman equations and ([7]) are compared with those our universe, expressed 
in the usual languages of the Einstein gravity with the "unknown" contributions of "dark 
energy" 3 , 



a 
a 



8ttG 
Ik 2 " 

AttG 
"Ik 2 " 



c 2 k 



(Pmatter + PD.E.) ~ > ' 

Kp 

matter + Pd.e.) + 3(p 

matter + PD.E.)]) 



(8) 
(9) 



2 In some literatures (see lOl lllj. for example), it has been claimed that, when considering perturbations 
around the background metric, the theory does not have an IR limit close to general relativity due to 
strongly coupled gravity fluctuations. This would be a very important issue for the consistency of the 
theory though the detailed discussions about this issue (see 12 1 for the troubles in [ll| ) is beyond the scope 
of the present work. But, in practice, this issue might not be quite relevant to our case since there would 
be a natural low momentum cut-off ~ ^/ 1 A^k I when considering fluctuations around cosmological solutions 
with a non- vanishing but tiny cosmological constant ~ h-w, as favored by the current observational data. 
This fact may change in our non-relativistic case since the meaning of the horizons would be also changed 
from the conventional ones, due to the momentum dependence of the light cones in UV. But in IR, i.e., low 
momentum, the usual meaning of the horizons would be "emerged" from the recovered Lorentz invariance 
(with A = 1) and so does the notion of the low momentum cut-off. Of course, we need some more rigorous 
analysis for a more explicit confirmation o f th is fact. 

3 I follow the physical convention of Ryden [lj| which disagrees with [l|, |2|: Ghoio = Gnorava/c 3 , Ah C ic = 

A r% 
-^Horava^ • 
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the energy density and pressure of the dark energy part can be read as ( for a related dis- 
cussion with matters in the context of the original Hofava gravity, see [151 ]) 

3k V / -k 2 2ku 2 \ 

PD - R " ± 8(3A-1) [rJo^ ~ W 2 ~ W ) ' (10) 

K 2 fi 2 ( k 2 Ikuj 2 \ 

PD - E - ~ T 8(3A^I) [b^ " W " "7 ' ( } 

respectively, where I have "defined" the fundamental constants of the speed of light c, the 
Newton's constant G, and the cosmological constant A as 

= «V|A W | G = A = 3 2 

8(3A-1) 2 ' 16tt(3A-1)' 2 w v ; 

Note that c 2 is non-negative always 4 , whereas G can be negative, i.e., anti-gravity, for 
A < 1/3, which implies that A c = 1/3 is the lower bound for the consistency with our 
universe. Moreover, with these definitions, one obtain the IR limit of the action ([2]) as the 
sum of the A-deformed Einstein-Hilbert action S'aeh and the IR limit of the dark energy 
action with 



1 /„ „„-„• . _o\ 2A 



Saeh = 87tG{ ^ x _ T) J dtd^N [1 (KyK* - \K 2 ) + - 



c 2 



(13) 



5d - K(IR) = 32^\| J dtd3x V9 NR{31 ( 14 ) 

The deformed action S'aeh agrees with the Einstein-Hilbert action when A = 1, but otherwise 
there are explicit A-dependances, generally 5 . This is in contrast to the directly measurable 
equations (jHJ), (JSJ) which define c and G, independently of A; actually, the Newton's constant 
in (jSJ) agrees with what we measure in laboratory, based on the Newton's Law of Gravity, 
F = —GMrn/R 2 14j and we can observe only the "renormalized" Newton's constant G, in 
contrast to [5J. 

Once the energy density and pressure of dark energy fflOl) . (fTT|) are identified, one can 
now compute the equation of state parameter as 

Pd.e. _ ( k 2 - 2kcua 2 - 3A 2 y a 4 \ 
WD - E - pd.e. \3k 2 + Qkuoa 2 + ) ' 1 j 



4 One might include oj term in the definition of c 2 , rather than including in pd.e., by replacing I A iy | —> 
\Aw\ i w but then c 2 can be negative when \uj\ > \A\y\. This is what has been considered in [5(] to get 
the non-vanishing speed of light for the flat universe limit — > with uj = ±8/i 2 (3A — 1)/k 2 , though 
our current universe seems to have a non- vanishing (positive) Aw But actually, there are infinitely many 
possible definitions of c 2 , depending on how much the u> term contributes to c 2 . In this paper, however, 
I do not consider all these possibilities and but consider only the simplest choice which is related to that 
of the Hofava's original proposal 0,13]. But this would be justified by experiments basically. 

and almost all other later works, S'aeh i s gi ven by S'aeh = 16 c 7rG / dtd 3 x^/gN[-^ {K^K^ — \K 2 ) + 
_ |A] with c 2 = re 4 Y A ^yl , G = A = ^A w c 2 , in the convention of this paper [jj]. This agrees 

with (fT3| for A = 1 but A < 1/3 is excluded from a mathematical consistency of c 2 > 0, rather than the 
physical reason of no anti-gravity, i.e., G > in the definitions of (TT!?|) . 
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FIG. 1: Plots of equation of state parameters wb.e. vs. scale factor a(t) for uj 2 < A^, kuo < 
( Q = -1/1.3,-1/2,-1/10, k = +1, or Q = +1/1.3, +1/2, +1/10, k = -1, with \A W \ = 1 (top 
to bottom in the left region )). When \uj\ is not far from |Avk|, there is a region where u?d.e. 
is fluctuating beyond the UV and IR limits. When \uj\ is small enough, wd.e. is monotonically 
decreasing from 1/3 in the UV limit to —1 in the IR limit. 



where I have introduced uj = ujRI, Aw = A^-Rq for the convenience. This interpolates from 
w^d.e. = 1/3 in the UV limit to w^.e. = ~ 1 in the IR limit but the detailed evolution patten 
in between them depends on the parameters k,uj,Aw- There are infinite discontinuities 
when uj 2 > A w , kuo < due to the vanishing Pd.e. with a non- vanishing pv.E. 0- But 
physically more interesting case would be uj 2 < A^?, kuj < or kuo > where there is no 
singular point of vanishing pd.e. but smoothly fluctuating/varying between the UV and IR 
limits (Fig. 1,2). For the original Hofava gravity with Uj = 0, yjp r is "always" monotonically 
decreasing from 1/3 in the UV limit to —1 in the IR limit [151 ]. 

In order to determine u>d.e., I need to know about the constant parameters k,uj,A w . 
Previously, those have been obtained by neglecting the matter contributions Q. Here, I 
consider the more improved analysis which does not need the consideration of the matters 
separately, based on the latest observational data. To this end, let me consider the series 
expansion of wd.e. in (fl5|) near the current epoch (a = 1), which coincides with Chevallier, 
Polarski, and Linder's parametrization exactly [161 ]. as 



wd.e. = w + w a (l - a) + w b (l - a) H (16) 



with 



k 2 - 2kuj - 3A 2 W _ 8k{ujk 2 + ujA 2 w + 2kA 2 w ) 
W ° ~ ?>{k 2 + 2kuj + A 2 w Y Wa ~ 3{k 2 + 2kuj + A 2 W ) 2 ' ( } 

So, by knowing wo and w a from the observational data, one can determine the constant 
parameters uj and A w as 

(1 — 2w — 3wq — w a )k 
(1 + 4w + Zwl + w a ) ' 
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FIG. 2: Plots of equation of state parameters u>d.e. vs. scale factor a(t) for kui > ( u = 
±2, ±l,±l/2 (top to bottom in the left region) for k = ±1, |A;y| = 1. In this case, u>d.e. is 
"always" monotonically decreasing from 1/3 in the UV limit to —1 in the IR limit. 



A 2 



(-1 + 9w 2 + 3w a )k 2 
3(1 + Awq + 3wi + w a ) 



(18) 



about the matter contributions separately, 
in which I have neglected the matter con- 
Kw from the transition point from de- 
the latest data sets (i.e., 
universe is allowed in the analy- 



from the 
clt in ©• From 
a non-flat 



Note that here I do not need to know 
in contrast to the previous analysis [~] 
tributions to get the approximate value of 
celeration phase to acceleration phase 
central values of the best fits) when 

ses 0, G1 ( w Q , w a ) = (-1.10,0.39), (-1.06,0.72), (-1.11,0.475), I get (ou,A w ) = 
(1.32,2.44), (1.14,2.10), (1.30,2.29), respectively and k = -1. Once the two constant 
parameters u,A w are determined, the whole function w D E .(a) is "completely" determined. 
From the obtained data, I plot the curves of wd,e.{ z ) v - s - the astronomer's variable of 
redshift z = 1/a — 1 in Fig. 3 and these correspond to those of Fig.l since \u\ < Aw- It 
is interesting to note that these curves appear to give similar results with a nearly model 
independent analysis of type la supernovae (SNe la) for < z < 0.6 and similar ten- 
dencies in other analyses of Gold data sets that use the parametrization ([TBI) even for higher 
redshifts (13] . 6 In addition, f|T8|) gives some constraints on the allowed values of wq and w a 
such that A w > for the consistency of our theory, w a > |(1 — 9wq), — 1 — 4w — 3w^ or 
w a < |(1 — 9Wq), — 1 — 4w — 3wq (See Fig. 4) . This seems to agree with observational 
data at about la (68.3%) confidence level 



20j, |2l|, |23j. This consistency condition may 



be tested in the near future, by sharpening the data sets. 

On the other hand, using the corresponding data f2& = —0.009, —0.000, 
in the current epoch (a = 1) for the deviation from the critical density, Qk = 



-0.0008 

1 



171,118 



6 In these analyses, the flat universe has been assumed but the results would be almost the same even for 
a non-flat universe. This is due to the fact that the relaxation of flatness broadens the ranges of wq and 



w a but the central values are almost unchanged 

0,0,0,0- 
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FIG. 3: Plots of equation of state parameters iod.e. vs - redshift z = 1/a — 1 for 
the latest data sets (u,K w ) = (1.32,2.44), (1.14,2.10), (1.30,2.29) from (w ,w o ) = 
(-1.10,0.39), (-1.06,0.72), (-1.11,0.475) (bottom to top) and k = — 1. 




H 2 k\Aw\Lp/2a 2 H 2 R 2 ] Mp, Hubble parameter H = a/a, the ratio of Planck mass and length 
Mp/Lp = c 2 /8nG, 1 get ^ — 0.0013,0.000,0.0004 (H R Mp/Lp) with the current value of 
Hubble parameter H . (See Table 1 for a summary of the data sets and their corresponding 
constant parameters, in the appropriate units.) 

In conclusion, I have considered the dark energy as a possible test of Hofava gravity. It 
seems that the IR modified Hofava gravity seems to be consistent with the current obser- 
vational data but we need some more refined data sets to see whether the theory is really 
consistent with our universe. However, it would be still useful to consider the list of possible 
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TABLE I: A summary of the data sets without assuming the flat universe in a priori and 
their corresponding constant parameters, in the conventional units of Hq (km s _1 Mpc _1 ) and 
M (H R M P /Lp). 



scenarios as follows. 

1. If k = 0, i.e., flat universe, as predicted by inflationary cosmology but not compulsory 
in the latest analyses (ItI 18 . 21 . 22| . is confirmed, there is no effect of the Hofava gravity 



in the FRW cosmology. But even in this case, its effect to the anisotropic cosmology and 
non-Gaussianity would be still open problems. 

2. If u>d.e. < — 1 an d k ^ 0, the original Hofava gravity with the detailed balance 
condition, which predicts —1 < u>d.e. < 1/3, may be ruled out. According to the current 
observational data, this scenario seems to be quite plausible and this is also consistent with 
other theoretical considerations [1, 0, d, 0] . 

3. Even if u>d.e. < ~~ 1 ; k ^ 0, and good agreements for small z are confirmed by 
determining the constant parameters to and Aw, some disagreements or inconsistencies for 
higher z can occur; for example, by comparing the higher-order term Wb in (jlfip with the 
experiments. In this case, one might consider several further modifications of the Hofava 
gravity by introducing other detailed-balance breaking terms with the additional constant 
parameters to control the disagreements. But one does not know how much new terms are 
needed minimally, in a priori 7 . Or, one might consider another definition of pd.e. and pb.e., 
by considering different definitions of the speed of light, rather than the simplest choice (fl~2l . 
as was discussed in the footnote no.l. 
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